Negative heat capacity at phase-separation in macroscopic systems. 
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Systems with long-range as well with short-range interactions should necessarily have a convex 
entropy S(E) at proper phase transitions of first order, i.e. when a separation of phases occurs. 

Here the microcanonical heat capacity c(E) = —^§Tgj^r is negative. This should be observable 
even in macroscopic systems when energy fluctuations with the surrounding world can be sufficiently 
suppressed. 



The argument is simple c.f . yj : At phase separation the 
weight e s ( E ' > ~ E / T of the configurations with energy E in 
the definition of the canonical partition sum 



Z(T) 



„S(E)-E/T dE 



becomes bimodal, at the transition temperature it has two 
peaks, the "liquid" and the "gas" configurations which 
are separated in energy by the latent heat. Consequently 
S(E) must be convex (d 2 S/dE 2 > 0, like y — x 2 ) and 
the weight in has a minimum at E m i n between the 
two pure phases. Of course, the minimum can only be 
seen in the microcanonical ensemble where the energy 
is controlled and its fluctuations forbidden. Otherwise, 
the system would fluctuate between the two pure phases 
by an, for macroscopic systems even macroscopic, energy 
AE ~ Ei a t oc N of the order of the latent heat in clear 
contrast to the usual assumption of the fluctuations in 
the canonical ensemble SE oc \f~N . The heat capacity is 



TABLE I: Parameters of the liquid-gas transition of small 
sodium clusters (MMMC-calculation 0, fl) i n comparison 
with the bulk for a rising number N of atoms, N sur f is 
the average number of surface atoms (estimated here as 



(1) J2 N ciuster) of a11 clusters with Ni > 2 together. a/T tr = 



As sur f*N/N sur f corresponds to the surface tension. Its bulk 
value is adjusted to agree with the experimental values of the 
a 3 parameter which we used in the liquid-drop formula for the 
binding energies of small clusters, c.f. Brechignac et al. 0, 
and which are used in this calculation [2| for the individual 
clusters. 





N 


200 


1000 


3000 


bulk 




T tr [K] 


940 


990 


1095 


1156 




qiat [eV] 


0.82 


0.91 


0.94 


0.923 


Na 


Sboil 


10.1 


10.7 


9.9 


9.267 




^<$surf 


0.55 


0.56 


0.44 






■^V surf 


39.94 


98.53 


186.6 


oo 




cr/Ttr 


2.75 


5.68 


7.07 


7.41 



X 



C v (E mm ) = 8E/dT = - 



(dS/dE) 2 
d 2 S/dE 2 



< 0. 



(2) 



I.e. the convexity of S(E) and the negative microcanoni- 
cal heat capacity are the generic and necessary signals of 
any phase- separation^. 

It is amusing that this fact, which is essential for the 
original purpose of Thermodynamics to describe boiling 
water in steam engines, seems never been really recog- 
nized in the past 150 years. However, such macroscopic 
energy fluctuations and the resulting negative specific 
heat are already early discussed in high-energy physics 
by Carlitz 0. 

This "convex intruder" in S(E) with the depth 
AS sur f(E) has a direct physical significance: Its depth 
is the surface entropy due to constraints by the existence 
of the inter-phase boundary between the droplets of the 
condensed phase and the gas phase and the correspond- 
ing correlation. AS sur f(E) is directly related to the sur- 



face tension per surface atom N sur f of the droplets. 

AS sur f (E m i n ^ 



/T tr 



N, 



surf 



(3) 



In my paper together with M.Madjet 3 we have com- 
pared the the values of AS sur f(E) calculated by Monte- 
Carlo with the use of a realistic interaction with the val- 
ues of the surface tension of the corresponding macro- 
scopic system. Table (UJshows the scaling behavior of 
AS sur f (E) with the size N of the system. 
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Roughly AS'su,,/ oc iV 2 / 3 and one may argue that this 
will vanish compared to the ordinary leading volume term 
Svoi(E) oc N. However, this is not so as S vo i(E) at en- 
ergies inside the phase-separation region (the convex in- 
truder) is the concave hull of S(E) (its slope gives the 
Maxwell construction of the caloric curve T(E)). It is a 
straight line and its curvature d 2 S vo i/dE 2 = 0. Conse- 
quently for large N 
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8 2 S/dE 2 ~ d 2 S voi /dE 2 + d 2 AS surf /dE 2 
x d 2 AS surf /dE 2 



(4) 



FIG. 1: Microcanonical Monte-Carlo (MMMC) @, Q simu- 
lation of the entropy s(e) per atom (e in eV per atom) of a 
system of N = 1000 sodium atoms at an external pressure 
of 1 atm. At the energy e < ei the system is in the pure 
liquid phase and at e > es in the pure gas phase, of course 
with fluctuations. The latent heat per atom is qi at =63 — 61. 
Attention: the curve s(e) is artificially sheared by subtracting 
a linear function 25 + e * 11.5 in order to make the convex 
intruder visible. s(e) is always a steep monotonic rising func- 
tion. We clearly see the global concave (downwards bending) 
nature of s(e) and its convex intruder. Its depth is the entropy 
loss due to additional correlations by the interfaces. It scales 
oc 7V -1//3 . From this one can calculate the surface tension 
per surface atom cr Bur f /Ttr = As sur f * N/N sur f. The double 
tangent (Gibbs construction) is the concave hull of s(e). Its 
derivative gives the Maxwell line in the caloric curve e(T) at 
Ttr- In the thermodynamic limit the intruder would disappear 
and s(e) would approach the double tangent from below, not 
of course S(E), which remains deeply convex: The probability 
of configurations with phase-separations is suppressed by the 
(infinitesimal small) factor e~ N 1 relative to the pure phases 
and the distribution remains strictly bimodal in the canonical 
ensemble. The region ei < e < e% of phase separation gets 
lost. 



and the depth of the intruder AS sur f = a/T tr * N sur f ~ 
TV 2 / 3 goes to infinity in the thermodynamic limit. Of 
course, the ubiquitous phenomena of phase separation 
exist only by this reason. 

It determines the (negative) heat capacity 



(dS/dEf 
d 2 S/dE 2 



< 0. 



(5) 



The physical (quite surprising) consequences are dis- 
cussed in P, |6j . 



Discussions with St.Ruffo are gratefully acknowledged. 
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